Scaling behavior in a quantum wire with scatterers 
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We study the conductance properties of a straight two-dimensional quantum wire with impurities 
modeled by s-like scatterers. Their presence can lead to strong inter-channel coupling. It was shown 
that such systems depend sensitively on the number of transverse modes included. Based on a poor 
man's scaling technique we include the effect of higher modes in a renormalized coupling constant. 
We therefore show that the low-energy behavior of the wire is dominated by only a few modes, which 
hence is a way to reduce the necessary computing power. The technique is successfully applied to 
the case of one and two s-like scatterers. 
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I. INTRODUCTION 

Over the past decade, quantum transport in two- 
dimensional electron systems has attracted increasing at- 
tention. In electron waveguide structures, the transmis- 
sion and resonance phenomena have been studied exten- 
sively, especially because the conductance was shown to 
be directly related to the scattering properties of the 
system. i-iTjhis relation is known as Landauer-Biittiker 
formulajjtj 



G =T T ' 



(1) 



where G denotes the 2-point conductance and T the full 
transmission function (spin degrees of freedom are ne- 
glected). Eq. (|l|) results in a quantized conductance for 
ballistic channels, which has been experimentally verified 
for a large range of conduction channels.^ 

It is therefore of great interest to study the properties 
of ballistic wires containing impurities. Especially for at- 
tractive scatterers the combined effect of the scatterer 
itself and the backscattering off the walls leads to inter- 
esting .-phenomena such as resonances and quasi-bound 
states.LTO Moreover strong inter-channel coupling is in- 
troduced into the system making calculations quite com- 
plicated. A simple model containing all of the features 
mentioned above consists of a pointJikp scatterer mod- 
eled by a two-dimensional (5-functionljcj In this paper we 
present a way to dramatically reduce the computational 
complexity of this problem by including the effect of the 
higher modes in a renormalized coupling constant. 

The outline of this paper is as follows: In the next 
section we present the standard approach to calculating 
the conductance in this kind of waveguide setup. We 



then briefly discuss results for two model systems for an 
impurity in the waveguide. For the case of a ^-function 
impurity, we introduce a scaling approach and discuss the 
results for the rescaled system. We finally generalize the 
scaling approach to a system with many scatterers and 
show results for the case of two <5-function scatterers. 



II. CONDUCTANCE CALCULATION 

In this section we sketch the calculation of the conduc- 
tance of an electron waveguide including a general scat- 
terer. We restrict ourselves to a single scatterer, because 
the transmission matrix for many scatterers can be ob- 
tained by multiplying together the transmission matrices 
for the individual single scatterers. 

Let us first obtain a solution of the Schrodinger equa- 
tion for an electron in a two-dimensional waveguide with 
a general scatterer. The Hamiltonian is given by 



H 



P 
2m 



+ V{x,y) + V c {y). 



(2) 



V c represents a confinement potential restricting the 
transverse movement of the electron and V(x,y) rep- 
resents the scattering potential. We assume V(x, y) to 
be non-zero only within a finite region small compared 
to the width of the channel. We can now expand any 
stationary solution y^E(x,y) of the Schrodinger equa- 
tion, HipE(x,y) — EipE(x,y), in a Fourier series with 
x-dependent expansion coefficients using the complete set 
of transversal modes, 



^E{x,y) 



^2c n (x) Xn (y), 



(3) 



where \n (y) is the transverse state of the electron in the 
absence of the scatterer. Inserting this series into the 
Schrodinger equation and employing orthogonality of the 
transversal modes we obtain a set of coupled equations, 

— -^c m (x) + k 2 m c m (x) = ^M mn (x)c n (x), (4) 

n 

with the definitions of the mode coupling constants 
M mn (x) and the wave vector k m , 



M mn (x) 



2m 



dyx* m (y)v(x,y)xn{y), 



2m 



k m — \j —j- (E — E n 



(5) 



(6) 



Outside the scattering region the solution is given by lin- 
ear combinations of the form 

, s _ J A n e iknX + B n e- lk " x x « , . 

Cn{X)- | Cnf ,ik„x + Dnf ,-ik„x x>0 K') 

For k n real, we get propagating modes, for k n imaginary 
evanescent modes. By matching the wavefunctions for 
both regions with the appropriate boundary conditions, 
the coefficients for incoming and outgoing waves can be 
obtained. 

The transmission coefficient for propagating modes is 



defined as T mn 
function as 



h I A I 1 



and the total transmission 



T(E) = Y, T " 



(8) 



where the sum extends over all propagating modes. The 
conductance is finally calculated using Eq. (|l|). 

The transmission coefficients can also be obtained from 
the Green's function of the system. More specifically, it is 
the retarded Green's function that governs this behavjp; 
and via its poles also describes resonance behavior .@'E£l 
The retarded Green's function for the wire without scat- 
terer is given by 



G° (x,y,x',y') = Y^n(y)Xn(y') 



,. 2m e 



ik n \x—x'\ 



n=l 



Z/hryif-i 



(9) 



or 



G ab (x, x') = Jdy dy> X ;(y)G° (x, y, x\ y 1 ) X b(y') (10) 

OC Sab- 

The full solution can then be obtained from the Dyson 
equation 

G a b\X 1 x ) = G ab {x, x ) 

+ ]T fdx" G° ac (x, x")V cd (x")G db (x", x') (11) 



with 



V cd (x")= dy X *MV(x",y)xd(y)- 



(12) 



III. SIMPLE MODEL SYSTEMS 



We will now apply the general outline of the previ- 
ous section to two simple model systems: the case of the 
(5-function scatterer described by 



V(x,y) = 7<S(x)% - yo) 



(13) 



and discussed by various authors JaldlLa as well as the 
slightly more realistic model with 



V(x,y) = 7<S(x)-^=-exp (-^ 



(14) 



describins-|a softer "impurity" scatterer in transversal 
direction.E3 

For both potentials, the waveguide can be split in two 
separate regions at x = and Eq. (0) can be written as 

( A P ik n x I f> -ik n x ^- n, 

C n [X) - < Cne ik nX + Dne ^k n x X>Q UOJ 



As ij) must be continuous at x = and its derivative must 
have a finite jump there, the same conditions must hold 
for the expansion coefficients c n (x). Thus using these two 
conditions on Eq. (Q) with the ansatz in Eq. (|l5|) yields 



A n + %- L^n + E) n , 

ik n (C„ — D n ) — ik n [A n — B n ) 

= Y M nm (A m + B rl 



(16) 

(17) 



If i\) is an evanescent mode, we can set k n — in n and 
must require A n — and D n = to have a normalizable 
wavefunction. The coupling constants are given by 



M, 



2m^i 



mil — 2 



x*„(yo)xm(yo) 



(18) 



for the case of the 5-function scatterer, Eq. (J13|), and 



M nm = — 2 -dyx* m (y)-^eiq>l ) Xn(y) (19) 



for the "impurity" scatterer, Eq. (|14|). In either case, the 
coupling constants do not depend on x anymore. For 
given coupling constants M nm , the transmission coeffi- 
cients can be computed by solving Eqs. (O) and dlj). 
For an attractive potential (7 < 0) resonance dips can be 
observed just below the energy at which a new propagat- 
ing mode develops. Their position and width depends 
sensitively-, on the number of modes and the scatterer 's 
strengthJia 



We note that in the case of the (5-function scatterer, the 
integral equation for the full Green's function, Eq. (11), 
can be solved explicitly to yieldEJ 



G(x,y,x',y') = G°(x,y,x',y') 

G o (x,y,0,y o )G o (0,yo,x',y') 
l/7-G°(0,yo,0,yo) 



(20) 



We first briefly summarize the results for the (5-function 
scatterer. As there is no solution to a two-dimensional 
(5-function the problem would be ill-posed unless one 
restricts the problem to a finite number of transverse 
modes. The potential is then equivalent to an s-like scat- 
terer. As more modes are included the electron becomes 
more strongly localized. For the following computations, 
we will use the parameters D — 300 A for the width of the 
channel (hard wall potential), yo — j^D for the transver- 
sal position of the scatterer, and the mass m = 0.067m e 
as the effective mass of an electron in GaAs-AlGaAs het- 
erostructures. The total number of modes will be de- 
noted by n c . For a potential strength of 7 = =F7 feVcm 2 , 
the result ing conductance is shown in Figs. VII (dashed 
line) and |VI1| (dotted line) . The dependence of the pole 
location on the number of modes can be seen in Fig. H 
(dots), with no indication of convergence with respect to 
n c . 

This convergence problem can be avoided by modifying 
the 5-function potential in y-direction to the smoother 
form of a Gaussian function as in Eq. (|l4|). The finite 
width allows the higher modes to decouple from the lower 
ones. In the limit p — > 0, we recover the (5-function poten- 
tial with the potential strength 7. By numerically eval- 
uating the integrals for the coupling constants, Eq. (19), 
one can again obtain the full transmission function and 
thus the conductance. As expected, the conductance 
curves closely resembles the results for the (5-function 
scatterer, showing the same characteristic resonance dip 
for attractive scatterers. Numerical results for the exact 
pole location of the transmission function are shown in 
Fig. [| using the potential parameters 7 = — 7feVcm 2 
and p = 3.0 A. In contrast to the (5-function potential, 
the result converges to its exact solution when using a 
large enough, but finite number of modes. 

Nevertheless, the characteristic behavior of the con- 
ductance is identical to the (5-function scatterer with a fi- 
nite number of modes. This again demonstrates that the 
5-function scatterer serves its purpose as a useful model 
if it is interpreted as an s-like scatterer for a fixed finite 
number of modes. 



IV. SCALING APPROACH 

The methods presented in the last section work fine 
except that they require an infinite or very large (for a 
strongly localized scatterer) number of modes to be in- 
cluded. Simply cutting of the sum at some small number 
does not give the correct answer. But calculations for 



large n c have the drawback that substantial computing 
power is needed. Therefore it would be useful if one could 
find a way to reduce the number of modes in the calcu- 
lation without neglecting their influence. A well known 
technique for this kind of problem is the poor man's scal- 
ing approach.EZl The idea is the following: Consider a 
typical relevant process which includes an intermediate 
virtual state. The sum over the virtual states shall have 
a cut-off. If the system exhibits scaling behavior one 
can infinitesimally reduce the cut-off and incorporate the 
change into a new effective coupling constant for a pro- 
cess not involving these virtual states. By repeating these 
steps one can integrate out a large portion of the states. 
For the present system this means that we investigate 
the influence of decreasing the mode cut-off n c on the 
coupling constant 7.E3 In the spirit of poor man's scaling 
in the Kondo problem we look at a diagram as shown 
in Fig. ra which represents a typical process in the wire. 
We will first discuss this in terms of a general scatterer 
before we show the detailed calculation for the (5-function 
scatterer from above. The contribution of the diagram is 
given by 

£ J dx"dx<" G° aa (x, x")V ac {x")G°M',x'") 

xV cb (x"')G° bb (x'",x'). (21) 



c=\ 



The cut-off n c stops the summation of the intermediate 
modes. Reducing the cut-off by one mode amounts to 
a difference that equals the n c 's summand. This change 
should be incorporated into the new coupling constant 
that is taken for a process with only one scattering event, 
which would have the following form 



dx"G° aa (x,x")V^(x")G° bb (x",x') 



(22) 



where the V^ b (x") denotes a potential with renormal- 
ized coupling constant. Assuming that although n c is 
clearly integer we can treat the difference as infinites- 
imal, and that the potential can always be written as 
Vab(x) = jV a b(x), we can formally write 



<h. 



(vaUxny 1 



x dx'"V anc (x")G» (x",x'")V ncb (x' )■ (23) 



This equation is in general difficult to solve, and couM 
potentially contain an explicit dependence on a and bxB 
We therefore show how to solve this equation for the <5- 
function scatterer with a finite number of modes as dis- 
cussed above. The contribution of the diagram is given 
by jG 1 ^ n (0, 0) 7. For large n c we can thus write the 
change in 7 as 



S') 

Sn r 



2 2 2m 2 



Tm c y Q 
D 



Z/hrZir, 



(24) 



For large n c we can approximate k n 
grate Eq. ( gj ) to obtain 



^ and inte- 



1 

1 



2m 

IK 
m 

l2~ 



-Ci 



27ry n c 



2" V D 

lnn r + const.. 



1 



■ In n, 



const. 



(25) 



where the last equality holds again for large n c and Ci(x) 
is the integral cosine function. Since we expect a small 
change of the cut-off not to change the result as a whole 
we can deduce the renormalization group equation 



/7i 2 7r\ _ (h 2 -K 

n c exp = n c cxp — — 

\ mj J \ 7717 



(26) 



With the help of Eq. (E6|) we can now calculate the con- 
ductance for a model with a large number of modes n c by 
including only n c modes and the renormalized coupling 



7 



H TT 



In S* 



(27) 



The cut-off h c must be chosen such that the approxi- 
mations made above are valid, i. e. the renormalization 
procedure will eventually break down for small n c . More- 
over, there is one strict physical criterion for a minimum 
lower bound to n c : To reproduce the resonance behav- 
ior of the full system, it is at least necessary to include. 
one evanescent mode in the renormalized calculation.Ej 
We would also like to point out that a divergence of 7 
is not only possible in Eq. (27), but furthermore that it 
may even abruptly change sign from —00 to +00. This is 
however not problematic: As can been seen from Eq. (20) 
it is I/7 that enters the Green's function, telling us that 
7 — ► ±00 show the same physics, independently of the 
sign of the coupling strength. With increasing scatterer 
strength, the transmittive behavior of the attractive scat- 
terer indeed nppre and more resembles that of the repul- 
sive scatterer B 



again plot the conductance for the attractive case. This 
time we scale down to n c = 50, 10, and 5. For large n c 
the method works very well, but as we get to small val- 
ues of n c , it eventually breaks down. However we used 
the assumption that n c is large in our derivation, so this 
breakdown is expected. But even for n c = 5 it still gives 
fairly good results. Another limit is reached when the 
number of propagating modes is close to n c . But only 
for n plop > n c it breaks down completely, as there is no 
evanescent mode left which could build up the resonance. 
It may seem surprising that such a simple method yields 
these results. However, it confirms that the type of pro- 
cess we have chosen is indeed the relevant one. 

An alternative description of resonance behavior can 
be given in terms of quasi-bound state. They are de- 
fined as the poles of the energy Green's function in the 
complex energy plane. Hence our approach should ap- 
ply to the quantities of this approach as well. In Fig. pi 
we show that the location of the poles can be calculated 
very well with the scaling method. It is remarkable that 
even down to n c — 2, all poles lie very close to the cor- 
rect value. Moreover, they are not randomly scattered 
around n c = 100, but rather lie on the parabolic curve of 
the poles for increasing n c . 

Our result shows that although a true description of 
quantum wires with scatterers requires the inclusion of 
all modes, the presence of higher modes does not have 
a strong affect on the low energy behavior. Their influ- 
ence can be cast in a renormalized value of the coupling 
constant. Therefore one can say that calculations done 
on these waveguide systems with a few number of modes 
should give qualitatively correct results. This is not a 
trivial result, because the mechanism of the suppression 
of the higher modes is not of thermal origin. Rather it 
is the scaling behavior of these systems that allows our 
conclusion. 



VI. MANY SCATTERER CASE 



V. COMPARISON TO EXACT FORMALISM 

In this section we show and compare the results when 
the formalism from above is applied to the single 6- 
function scatterer given by Eq. (fL3|). At the end we dis- 



cuss the implications of this scaling behavior. In Fig. VII 
we show the conductance of a system with an attractive 
(top) and a repulsive (bottom) (^-function scatterer with 
a finite number of modes. The shape of the curves is 
altered by a variation of n Cl especially the width and po- 
sition of the resonance in the attractive case are quite 
sensitive (comparing e. g. the dotted and dashed lines). 
When we apply the scaling approach, Eq. (E71) to the 
model with n c — 100 and scale it down to n c = 10, we 
obtain an almost identical result (solid line) . This clea rly 
confirms the validity of our approach. In Fig. VH we 



Another important consequence is that future calcu- 
lations need to be done only for a few low lying modes. 
This is important especially for systems with many scat- 
terers: Due to the scattering matrix multiplication, the 
computing time grows linearly with the number of scat- 
terers and grows more than quadratic with the number 
of modes. When choosing n c one needs of course to make 
sure, that the evanescent modes are strongly enough lp=. 
calized, such that they do not carry a virtual currentE2l 
and thereby provide a minimum upper cut-off. In order 
to demonstrate this, we perform a calculation for the sim- 
ple case of two (^-function scatterers placed in series and 
separated by a distance d. For a few number of modes 
this model has been discussed in Reference |2l]. Here we 
present results for the case of many modes and discuss 
how they can be obtained using our scaling method. 

The specific system under consideration consists of two 



S- function scatterers located at (0, y\) and (d,y 2 ), with 
potential strengths 71 and 72 respectively. The waveg- 
uide itself shall be the same as before. The Green's func- 
tion of this system can be obtained analytically and is 
given by the lengthy however very useful expression 



G(x, y, x', y') = G°(x, y, x', y') + 7l G (x, y, 0, yi ) 

+ j 2 G°{x,y,d,y 2 ) 



1 — aia-2 
where the following abbreviations have been used 



c\ + a\c 2 
1 - aia 2 

(28) 



C\ 



(>2 



«1 



«2 



G°(0, yi ,x>, y >) 
l-TiGOfO.ift.O.i/i)' 

G°(d : y 2 ,x>,y>) 
l- l2 G°(d,y 2 ,d,y 2 ) : 

l2 G°(0, yi ,d,y 2 ) 
l- 7 iG (0,2/i,0, yi ) ; 

7ig°(d,y 2 ,0,yi) 
l-7iG0(d,j/ 2 ,d,t/ 2 )' 



(29) 



To apply our scaling method, we employ the approxi- 
mation that the evanescent modes fall off exponentially. 
Hence contributions to the scaling of order 7172 are 
strongly suppressed for higher modes, i. e. |fc n d| is large. 
We can then simply use the single scatterer result and 
apply it independently t o bo th coupling constants. The 
results ar e sho wn in Fig. VII a) and b). 

In Fig. VII a) we compare the conductance for vari- 
ous number of modes. Already for two modes the peaks 
in the lowest subband with one evanescent mode are 
present. The inclusion of higher modes however leads to 
a strong shift in position and shape. Although the true 
peak structure becomes visible for n c — 10, it is fully 
developed only for many more modes. In addition, the 
results in the second subband are completely wrong for 
two modes. This should be ex pected since no evanescent 
mode is present. In Fig. VII b) we compare our scaling 
result to the exact one and again find that it works very 
well. It is remarkable that already for two modes the 
deviation in the first subband is negligible. 

The generalization to many scatterers is obvious: one 
only has to keep in mind to take enough evanescent 
modes to mimic the decaying currents between the scat- 
terers. The Fabry-Perot behavior that causes the oscil- 
lations is created by the propagating modes, however it 
can be enhanced by the first few evanescent ones. 



VII. CONCLUSION 

In this article we discussed the conductance proper- 
ties of ballistic quantum wires containing impurities. As 
shown, they depend strongly on the number of modes 
included and the calculations can thus become quite dif- 
ficult. We proved that it is possible to incorporate the 



effect of higher modes into an effective coupling constant 
for the lower lying modes. This was done by a poor man's 
scaling approach. For the (5-function scatterer with a 
fixed number of modes we compared our result to exact 
calculations and found very good agreement. We there- 
fore conclude that our approach does not only allow for 
quicker calculations but also shows that the system's be- 
havior is governed by a few low-lying modes. 
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FIG. 3. Poles of the Green's function as a function of n c . 
The dots indicate the pole location for values from n c = 2 to 
n c = 110. The diamonds represent the case n c = 100 with n c 
varying between 2 and 60. 
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FIG. 1. The conductance for an attractive (top) and a re- 
pulsive (bottom) 5-scatterer of strength 7 = =p7 feV cm 2 and 
n c = 100 (dashed), n c = 10 (dotted) and n c = 100 scaled 
down to n c = 10 (solid). The original curve with n c = 100 
and the rescaled one with n c = 10 lie nearly on top of each 
other. 
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FIG. 4. Poles of the Green's function as a function of n c . 
The circles indicate the pole location for values from n c = 2 
to lie — 110. When using n c = 80 modes, the pole location is 
converged to its final position. 
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FIG. 5. A diagram describing the propagation from a state 

characterized by mode n and longitudinal momentum k via 

two scattering events into the state n',k' . 
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FIG. 2. The conductance for an attractive (5-scatterer of 
strength 7 = — 7feVcm 2 , n c = 100 calculated without scal- 
ing (dotted) and with n c rescaled to 10 (solid) , to 50 (dashed) , 
and to 5 (dot-dashed). If the number of modes after the scal- 
ing is large enough, the difference is hardly visible. 
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FIG. 6. Conductance for the 2 <5-function scatterer 
case. The used parameter values are 71 = 25feVcm , 
72 = —25 feV cm 2 and d = 90 nm. In a) the different lines are 
n c — 50 (solid), 30 (dotted), 10 (dashed) and 2 (dot-dashed). 
One can see that already from n c = 10 on, the qualitative 
behavior is the same. 

In b) the scaling is applied to the n c = 100 case (solid) 
and performed down to n c = 10 (dotted) and 2 modes 
(dot-dashed) . The dashed line shows how well it compares 
to the case with n c = 10 modes without scaling. 



